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Abstract 

Inclusive jet cross sections for events with a large rapidity gap with respect to the 
proton direction from the reaction ep — > jet + X with quasi-real photons have been 
measured with the ZEUS detector. The cross sections refer to jets with transverse energies 
E 3 ^ 1 > 8 GeV. The data show the characteristics of a diffractive process mediated by 
pomeron exchange. Assuming that the events are due to the exchange of a pomeron 
with partonic structure, the quark and gluon content of the pomeron is probed at a scale 
~ (E^ 1 ) 2 . A comparison of the measurements with model predictions based on QCD plus 
Regge phenomenology requires a contribution of partons with a hard momentum density in 
the pomeron. A combined analysis of the jet cross sections and recent ZEUS measurements 
of the diffractive structure function in deep inelastic scattering gives the first experimental 
evidence for the gluon content of the pomeron in diffractive hard scattering processes. The 
data indicate that between 30% and 80% of the momentum of the pomeron carried by 
partons is due to hard gluons. 
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1 Introduction 



Electron-proton collisions at HERA have shown evidence for hard processes in diffractive re- 
actions. Both in deep inelastic scattering (DIS) {Q 2 > 10 GeV 2 , where Q 2 is the virtuality of 
the exchanged photon) |], |], Q and in photoproduction (Q 2 ~ 0) j|, |5j, events characterized 
by a large rapidity gap towards the proton direction have been observed and interpreted as 
resulting from diffractive scattering M, |2|, |J. In the DIS regime, hard scattering for this class 
of events has been revealed through the virtuality of the probing photon [Q, [J and through 
the observation of jet structure in the final state 0. In the photoproduction domain, the hard 
scattering has been identified through jet production [|, [|. 

Diffractive processes are generally considered to proceed through the exchange of a colour- 
less object with the quantum numbers of the vacuum, generically called the pomeron (iP). 
Although the description of soft diffractive processes in terms of pomeron exchange has been 
a phenomenological success, the description of the pomeron in terms of a parton structure at 
first lacked experimental support. On the basis of pp data || from the CERN ISR, Ingelman 
and Schlein |7j] suggested that the pomeron may have a partonic structure. The observation of 
jet production in pp collisions with a tagged proton (or antiproton) made by the UA8 Collab- 
oration H gave strong evidence for such a structure. Further evidence has been provided by 
the observations made at HERA [1 — 5], which in addition include the first measurements of 



the diffractive structure function in DIS |9|, [Hj 



The description of diffractive processes in terms of QCD has remained elusive, in part due 
to the lack of a sufficiently large momentum transfer on which to base the perturbative expan- 
sion. Cross sections for diffractive processes involving large transverse energy jets or leptons 
in the final state are, however, amenable to perturbative QCD calculations [7, 11 — 17]. Their 
measurement could answer several questions concerning the structure of the pomeron such as: 
whether the parton picture is valid for the pomeron and universal parton densities can be de- 
fined; what fraction of the pomeron momentum is carried by gluons and what by quarks; and 
whether a momentum sum rule applies to the pomeron 



This paper presents the first measurement of inclusive jet cross sections in photoproduction 
at centre-of-mass energies ~ 200 GeV with a large rapidity gap. This process is sensitive to both 
the gluon and quark content of the pomeron. In order to examine the partonic structure of the 
pomeron, these jet cross sections are compared to predictions from models based on perturbative 
QCD and Regge phenomenology. The jet cross sections measured in photoproduction, combined 
with the results on the diffractive structure function in deep inelastic scattering [Kj, give the 



first experimental evidence for the gluon content of the pomeron. The result does not depend 
on the flux of pomerons from the proton nor on the assumption that a momentum sum rule can 
be defined for the pomeron. The data sample used in this analysis corresponds to an integrated 
luminosity of 0.55 pb _1 and was collected during 1993 with the ZEUS detector at HERA. 



2 Experimental setup 
2.1 HERA operation 

The experiment was performed at the electron-proton collider HERA using the ZEUS detector. 
During 1993 HERA operated with electrons of energy E e = 26.7 GeV colliding with protons of 
energy E p = 820 GeV. HERA is designed to run with 210 bunches separated by 96 ns in each of 
the electron and proton rings. For the 1993 data-taking, 84 bunches were filled for each beam 
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and an additional 10 electron and 6 proton bunches were left unpaired for background studies. 
The electron and proton beam currents were typically 10 mA, with instantaneous luminosities 
of approximately 6 • 10 29 cm -2 s _1 . 



2.2 The ZEUS detector and trigger conditions 

ZEUS is a multipurpose magnetic detector. The configuration for the 1993 running period 
has been described elsewhere pi |l9]j . A brief description concentrating on those parts of the 
detector relevant to this analysis is presented here. 

Charged particles are tracked by two concentric cylindrical drift chambers, the vertex de- 
tector (VXD) and the central tracking detector (CTD), operating in a magnetic field of 1.43 T 
provided by a thin superconducting coil. The coil is surrounded by a high-resolution uranium- 
scintillator calorimeter (CAL) divided into three parts, forward[| (FCAL) covering the pseudora- 
pidity^ regionP] 4.3 > rj^ > 1.1, barrel (BCAL) covering the central region 1.1 > 77^ > —0.75 
and rear (RCAL) covering the backward region —0.75 > rj^ > —3.8. The solid angle cover- 
age is 99.7% of 47T. The CAL parts are subdivided into towers which in turn are subdivided 
longitudinally into electromagnetic (EMC) and hadronic (HAC) sections. The sections are sub- 
divided into cells, each viewed by two photomultiplier tubes. The CAL is compensating, with 
equal response to hadrons and electrons. Measurements under test beam conditions show that 
the energy resolution is ue/E = 0.18/v^ (E in GeV) for electrons and <7e/E = 0.35/VE for 
hadrons [pOR . In the analysis presented here, CAL cells with EMC (HAC) energy below 60 MeV 
(110 MeV) are excluded to minimize the effect of calorimeter noise. This noise is dominated by 
the uranium activity and has an r.m.s. value below 19 MeV for EMC cells and below 30 MeV 
for HAC cells. For measuring the luminosity as well as for tagging very small Q 2 processes, two 
lead-scintillator calorimeters | 2l"| , located at 107 m and 35 m downstream from the interaction 



point in the electron direction, detect the bremsstrahlung photons and the scattered electrons 
respectively. 

Data were collected using a three-level trigger [JT9] . The first-level trigger (FLT) is built as 



a deadtime-free pipeline. The FLT for the sample of events analysed in this paper required 
a logical OR of different conditions on sums of energy in the CAL cells. The average FLT 
acceptance for the events under study was approximately 90%. The second-level trigger used 
information from a subset of detector components to differentiate physics events from back- 
grounds consisting mostly of proton beam gas interactions. The third-level trigger (TLT) used 
the full event information to apply specific physics selections. For this analysis, the following 
conditions were required: a) the event has a vertex reconstructed by the tracking chambers 
(VXD+CTD) with the Z value in the range \Z\ < 75 cm; b) E — pz > 8 GeV, where E is the 
total energy as measured by the CAL 



pz is the Z-component of the vector 



£ = ££*, 



1 The ZEUS coordinate system is defined as right-handed with the Z axis pointing in the proton beam 
direction, hereafter referred to as forward, and the X axis horizontal, pointing towards the centre of HERA. 

2 The pseudorapidity is defined as — ln(tan|), where the polar angle 9 is taken with respect to the proton 
beam direction, and is denoted by r)d (77) when the polar angle is measured with respect to the nominal interaction 
point (the reconstructed vertex of the interaction) . 

3 The FCAL has the forward edge at ?7edge = 4.3 with full acceptance for r/d < 3.7. 
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the sums run over all CAL cells, Ei is the energy of the calorimeter cell i and r*j is a unit 
vector along the line joining the reconstructed vertex and the geometric centre of the cell i; c) 
Pz/E < 0.94 to reject beam-gas interactions; and d) the total transverse energy as measured 
by the CAL, excluding the cells whose polar angles are below 10°, exceeds 12 GeV. 

3 Diffractive hard photoproduction 

Diffractive hard photoproduction processes in ep collisions are characterized by Q 2 ~ and 
by a final state consisting of a hadronic system X containing one or more jets, the scattered 
electron and the scattered proton 



where {pf) denotes the initial (final) state proton and X consists of at least one jet plus the 
remaining hadronic system (X r ). 

The kinematics of this process are described in terms of four variables. Two of them describe 
the electron-photon vertex and can be taken to be the virtuality of the exchanged photon (Q 2 ) 
and the inelasticity variable y defined by 



where E' e denotes the scattered electron energy and 9' e is the electron scattering angle. The 
other two variables describe the proton vertex: the fraction of the momentum of the initial 
proton carried by the scattered proton (#/), and the square of the momentum transfer (t) 
between the initial and final state proton. In terms of these variables and at low values of Q 2 
and t, the square of the mass of the hadronic system X is given by 



where s is the square of the ep centre-of-mass energy. 

Diffractive processes in which the photon dissociates give rise to a large rapidity gap between 
the hadronic system X and the scattered proton: 



going hadron belonging to the system X. The same signature is expected for double dissociation 
where the scattered proton is replaced by a low mass baryonic system (N). In this paper, the 
outgoing proton (or system N) was not observed, and instead of y^ x the pseudorapidity (i]^ x ) 
of the most forward-going hadron in the detector was used. 

Two cross sections are presented in this paper. First, the cross section for inclusive jet 
production is measured as a function of the pseudorapidity of the jet (rf et ) (for the definition of 
the jet variables see section 4) in reaction (1) with the most-forward going hadron at i]^ x < 1.8. 
This corresponds to a rapidity gap of at least 2.5 units measured from the edge of the CAL 
(Ar] GAP = r] cdgc — r\ h r nax)- This cross section is denoted by 



e + Pi — > e + X + pf — > e + (jet + X r ) + p f 



(1) 



E' e 1 - cos 6' e 
E~ e 2 



M\ « (1 - x f ) y s 



(2) 




(4) 



and is measured in the r] 3et range between —1 and 1. Second, the integrated cross section for 
inclusive jet production is determined as a function of rf max 

h d f +1 ' t d& h d 

a (Vmax < Vmax) = J ^ dlf e — — (r) max < T] max ) (5) 

and is measured in the range of f] max between 1 and 2.4. Both measurements include contribu- 
tions from double dissociation where the large-rapidity-gap requirement is satisfied. 

The jet cross sections refer to jets at the hadron level with a cone radius, R = ^/Ar] 2 + A0 2 , 
of one unit in pseudorapidity (rj) — azimuth ((f)) space and integrated over the transverse 
energy of the jet E 3 ^ 1 > 8 GeV. They are given in the kinematic region Q 2 < 4 GeV 2 and 
0.2 < y < 0.85. This region corresponds to photoproduction interactions at centre-of-mass 
energies in the range 130-270 GeV with a median Q 2 w 10~ 3 GeV 2 . 



3.1 Models 

The description of diffractive hard processes in terms of QCD is still in an early stage. Two 
main theoretical approaches have been considered. Both assume that a pomeron (IP) is emitted 
by the proton. The variable xjp = l— x/is then the fraction of the initial proton's momen- 
tum carried by the pomeron and M\ »s x^ys is the square of the 7-/P centre-of-mass energy. 
The two approaches differ in the modelling of jet production in 'jIP collisions. One of them 
[7, 11 — 13] assumes factorisation (factorisable models) while the other one does not [14, 16] 
(non-factorisable models). The latter are, however, not considered in what follows due to the 
lack of a Monte Carlo generator with an appropriate description of the event jet structure. 

Calculations based on factorisable models involve three basic ingredients: the flux of pomerons 
from the proton as a function of xjp and t, the parton densities in the pomeron and the matrix 
elements for jet production. The pomeron is assumed to be a source of partons which inter- 
act either with the photon (direct component) or with a partonic constituent of the photon 
(resolved component). As an example, the contribution of the direct component to the cross 
section for reaction (1) is given by 



Vdir = J dyf j/e (y) J J dx F dtfjp/ p (x]p,t)Y^ J d Pj2 J d PT — ^pP^ ( g » Pr» A* 2 ) fi/s>(P>t* 2 ) 



(6) 

where f 1 / e is the flux of photons from the electron^ and fjp/ p is the flux of pomerons from 
the proton. The sum in i runs over all possible types of partons present in the pomeron, and 
fi/jp(/3, /i 2 ) is the density of partons of type i carrying a fraction (3 of the pomeron momentum at 
a scale [i 2 and is assumed to be independent of t. The sum in j and k runs over all possible types 
of final state partons and a,i +1 ^j + k is the cross section for the two-body collision % + 7 — >• j + k 
and depends on the square of the centre-of-mass energy (s), the transverse momentum of the 
two outgoing partons (pr) and the momentum scale (//) at which the strong coupling constant 
(a s (fi 2 )) is evaluated. One possible choice is /i 2 = p\. In these models, the pomeron flux 
factor is extracted from hadron-hadron collisions using Regge theory, and the matrix elements 
are computed in perturbative QCD. However, the parton densities in the pomeron have to 
be extracted from experiment. While the recent measurements of the diffractive structure 
function in DIS at HERA |g [10 give information on the quark densities in the pomeron, the 



gluon content has so far not been determined. 

4 The Q 2 dependence has been integrated out using the Weizsacker- Williams approximation. 
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Two forms of the pomeron flux factor are commonly used. The Ingelman-Schlein form (IS) 
22 uses a parametrisation of UA4 data [2~3[|: 



fip/pixjp, i) = • (3.19e 8t + 0.212e 3t ) (7) 

Xjp 

where cq = ^3 GeV -2 and t is in GeV 2 . The Donnachie-Landshoff form (DL) |12| is calculated 
in Regge theory, with parameters determined by fits to hadron-hadron data: 

M P (^,t) = 0F 1 (t) 2 4- 2 ^ (8) 

using the elastic form factor Fi(t) of the proton, the pomeron-quark coupling 6 — 1-8 GeV -1 
and the pomeron trajectory a(t) = 1.085 + 0.25t with t in GeV 2 . 

Various parametrisations of the parton densities in the pomeron have been suggested on 
theoretical grounds [7, 11 — 13]. The following represent extreme possibilities for the shape of 
the quark and gluon momentum densities: 

• hard gluon density /3f g /ip(l3,fi 2 ) = 6/5(1 — (3); 

• soft gluon density /3f g /jp((3, /z 2 ) = 6(1 — /3) 5 ; 

• hard quark density (for two flavours) Pf q /]p(P, /i 2 ) = |/5(1 — 13). 

The first two assume a pomeron made entirely of gluons and the last one a pomeron made of 
uu and dd pairs. In all cases a possible /i 2 dependence of the parton densities is neglectedF] and 
the densities are normalised such that all of the pomeron's momentum is carried by the partons 
under consideration, £jp(/z 2 ) = J d(3J2i Pfi/jp(@, A* 2 ) — 1- 

However, since the pomeron is not a particle, it is unclear whether or not the normalisations of 
the pomeron flux factor and the momentum sum of the pomeron can be defined independently. 
Nevertheless, for an assumed normalisation of the flux factor, the momentum sum Sjp(// 2 ) can 
be measured. The definition used for the DL form of the pomeron flux factor is the appropriate 
one if the pomeron were an ordinary hadron and, hence, the one in which the momentum sum 
rule might be fulfilled [^4|]. 

Factorisable models presently account only for diffractive hard processes in which the proton 
remains intact. Since the measurements are based on the requirement of a large rapidity gap in 
the central detector, the contribution to the measured cross sections from double dissociation 
has to be taken into account when comparing with model predictions. 



4 Data selection and jet search 

Events from quasi-real photon proton collisions were selected using the same criteria as reported 
earlier ||25| . The main steps are briefly discussed here. 

Events satisfying the TLT selection described in section 2.2 are first selected. A cone algo- 



rithm in rj-(f) space with a cone radius of 1 unit |26|, [27]] is then used to reconstruct jets, for both 
data and simulated events (see next section) from the energy deposits in the CAL cells (cal 
jets), and for simulated events also from the final state hadrons (had jets). The axis of the jet is 

5 The /i 2 dependence of the parton densities in the pomeron is expected to be smaller than the differences 
between the various parton densities considered . 
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defined according to the Snowmass convention [ 57| , i] jet (0 je *) is the transverse energy weighted 
mean pseudorapidity (azimuth) of all the objects (CAL cells or final state hadrons) belonging 



to that jet. This procedure is explained in detail elsewhere | 25| . The variables associated with 
the cal jets are denoted by E^ cal , rj^, and while the ones for the had jets by ESjf , r) jet , 
and je *. 

A search for jet structure using the CAL cells is performed in the data. Events with at least 
one jet fulfilling the conditions E^^ > 6 GeV and — 1 < rf^ < 2 are retained. Beam-gas 
interactions, cosmic-ray showers, halo muons and DIS neutral current events are removed from 
the sample as described previously . The sample thus consists of events from ep interactions 
with Q 2 < 4 GeV 2 and a median Q 2 « 10~ 3 GeV 2 . The 7p centre-of-mass energy (W) is 
calculated using the expression W = ^Jys. The event sample is restricted to the kinematic 
range 0.2 < y < 0.85 using the following procedure. The method of Jacquet-Blondel [28| is 
used to estimate y from the energies measured in the CAL cells (see section 2.2) 

E-pz 

vjb - 

As can be verified using photoproduction events with an electron detected in the luminosity 
monitor (tagged events), yjs systematically underestimates y by approximately 20%, which 
is adequately reproduced in the Monte Carlo simulation of the detector. To allow for this 
effect, the event selection required 0.16 < yjs < 0.7. The sample thus obtained consists of 
19,485 events containing 24,504 jets. The only significant background, which is from misiden- 
tified DIS neutral current interactions with Q 2 > 4 GeV 2 , is estimated to be below 2%. The 
photoproduction origin of the sample is verified by the expected contribution (26%) of tagged 
events. 



5 Monte Carlo simulation 



Events from diffractive hard photoproduction processes were simulated using the program 
POMPYT [22|. These events were used to determine both the response of the detector to 
the hadronic final state and the correction factors for the cross sections for jet production with 
a large rapidity gap. 

The POMPYT generator is a Monte Carlo implementation of the model proposed in 0. 
The generator makes use of the program PYTHIA |29|] to simulate electron-pomeron inter- 
actions via resolved and direct photon processes. In PYTHIA, the lepton-photon vertex is 
modelled according to the Weizsacker- Williams approximation and the effects of initial state 
bremsstrahlung from the electron are simulated by using the next-to- leading order electron 
structure function [30]. Radiative corrections in our kinematic region, where W is larger than 



100 GeV, are expected to be negligible |JI|]. For the resolved processes, the parton densities of 



the photon were parametrised according to GS-HO |32] and evaluated at the momentum scale 
set by the transverse momentum of the two outgoing partons, fi 2 = p\. The parton densities 
in the pomeron were parametrised according to the forms described in section 3.1 and were 
taken to be independent of any scale. In PYTHIA, the partonic processes are simulated us- 
ing leading order matrix elements, with the inclusion of initial and final state parton showers. 



Fragmentation into hadrons was performed using the Lund string model |33[ as implemented 
in JETSET p3 . Samples of events were generated with different values of the minimum cutoff 
for the transverse momentum of the two outgoing partons, starting at pTmin — 3 GeV. 
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The program PYTHIA was also used to simulate standard (non-diffr active) hard photo- 
production events via resolved and direct photon processes. Events were generated using the 
leading order predictions of GRV for the photon parton densities and MRSD_ [|36| for the 
proton parton densities. 

All generated events were passed through the ZEUS detector and trigger simulation pro- 
grams. They were reconstructed using the same standard ZEUS off-line programs as for the 
data. 



6 Event characteristics 



The event variable rj max , as in previous studies by ZEUS [|], 0, £|, [HJ, was used to select events 
with a large rapidity gap. For the data, this variable is defined as the pseudorapidity (r]™ l ax ) 
of the most forward condensate with an energy above 400 MeV. A condensate is a contiguous 
energy deposit above 100 MeV for pure EMC and 200 MeV for HAC or mixed energy deposits 
in CAL. In the samples of simulated events, the f] m ax variable is defined at both the hadron 
and CAL levels. At the hadron level, all particles with lifetimes larger than 10 -13 s, energies 
in excess of 400 MeV and pseudorapidities below 4.5 are considered as candidates for the most 
forward final state particle, and T]^f x defines the pseudorapidity of the most forward particle. 
The CAL level uses the same definition as for the data. 

The mass of the hadronic system (Mx) of each event is reconstructed using the CAL cells, 
Mx l = \/E' 2 — p 2 The correlation between M™ 1 and rj™ l ax for the sample of events with 
at least one cal jet fulfilling the conditions E 3 ^ l cal > 6 GeV and — 1 < rfj% < 1 is displayed in 
Fig. |l]a. As shown in our previous publication [|J, there exists a distinct class of events with low 
Vmax values. The large-rapidity-gap events ijf^ ax < 1.8) are found to populate the region of low 
M™ 1 values, in contrast to the bulk of the data which have large M™ 1 values. These features 
of the data are reproduced by the Monte Carlo simulations: the events from a simulation of 
standard hard photoproduction processes using PYTHIA populate the region of large rj^ix an d 
large Mjjf' values; the events from a simulation of diffractive hard processes using POMPYT 
extend into the region of low r)^ l ax and Mjf ' values. 

A study of the region of low M™ 1 in the data sample reveals the following features. The 
Vmax distribution for events with M^ 1 < 30 GeV is shown in Fig. [TJb along with the predictions 
of PYTHIA and of POMPYT with a pomeron made of hard gluons (normalised to the number 
of data events above and below rf^ ax = 2.5, respectively). The simulation of non-diffractive 
processes by PYTHIA cannot reproduce the shape of the measured v^ ax distribution. On the 
other hand, the predictions of POMPYT describe well the shape of the data below if^ ax ~ 3. 

The Mx 1 distribution for the sample of events with if^ ax < 1.8 is shown in Fig. |I|c. This 
sample consists of 49 events containing 68 jets. The data exhibit an enhancement at low masses, 
15 GeV~ Mx 1 ~ 30 GeV, which is reproduced by the simulation of POMPYT with a hard gluon 
density (normalised to the number of data events). The W of each event is reconstructed using 
Djbi W cal = ^Ujb s - The distribution of W cal for events with v^ ax < 1.8 is shown in Fig. |T]d 
along with the expectations of POMPYT with a hard gluon density (normalised to the number 
of data events). The W cal dependence exhibited by the data sample is well reproduced by the 
simulation of POMPYT. The expectations of POMPYT using a hard quark density (not shown) 
also give a good description of the distributions of the data. Note that POMPYT assumes the 
cross section to be independent of W as expected for diffractive processes mediated by pomeron 
exchange. The good agreement with the data gives evidence for the diffractive nature of the 
large-rapidity-gap events. 
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In summary, the data exhibit a different behaviour in the region of low masses of the hadronic 
system compared to that of high masses. At low masses, the shape of the rf^ ax distribution 
in the data sample cannot be accounted for by the simulation of non-diffractive processes as 
in PYTHIA. On the other hand, the features of the data are described by the predictions 
of diffractive processes mediated by pomeron exchange as in POMPYT. These facts support 
the interpretation of these large-rapidity-gap events as being produced by diffractive processes 
via pomeron exchange. Therefore, the measurements of jet cross sections presented in the 
next section are compared to the predictions of models based on pomeron exchange. However, 
without a detected fast proton in the forward direction, the jet cross sections refer to events 
with a large rapidity gap. These include events with a diffractively scattered proton as well as 



those with a diffractively dissociated proton with mass less than approximately 4 GeV ||10|| . In 
this way, the measurements are presented in a model independent form suitable for comparison 
with calculations other than those presented here. 

6.1 Energy and acceptance corrections 

The method to correct the transverse energy of a jet as reconstructed using the CAL cells has 
been discussed elsewhere [EIJ. For samples of simulated events, the transverse energy of a jet as 
measured by the CAL (E^ cal ) was compared to that reconstructed using the final state hadrons 
(E^ )■ The corrections to the jet transverse energy were constructed as multiplicative factors, 
C '(E^aii rfcai) i which, when applied to the E T of the cal jets, give the corrected transverse 
energies of the jets: E^ 1 = C(E^ t cah r] :> c f l ) x E^lai- The function C corrects for energy losses, 
and for values E^lai > 10 GeV is approximately flat as a function of -E^Li an d varies between 
1.08 and 1.18 depending on rf^. For E^ l cal near threshold, E^ l cal m 6 GeV, this correction 
procedure can give values as large as 1.40. No correction is needed for rf et {rf et m rf^). The 
procedure was validated by comparing the momenta of the tracks in the cal jet in data and in 
Monte Carlo simulations. From this comparison it was concluded that the energy scale of the 



jets is corrected to within ± 5% p5| . The correction procedure was applied to the data sample 
of jets with E^ l cal > 6 GeV to select for further study those jets with corrected transverse 
energies of E^ 1 > 8 GeV and with the jet pseudorapidity in the range between —1 and 1. 

The events generated by POMPYT were used to compute the acceptance correction for 
the inclusive jet distributions. This correction function takes into account the efficiency of 
the trigger, the selection criteria and the purity and efficiency of the jet and r]^ x selection. 
It also corrects for the migrations in the variable rj^ l ax and yields cross sections for the true 
rapidity gap defined by i]^ x and rj = 4.5. After applying the jet transverse energy corrections, 
the purity was ~ 40% and the efficiency was ~ 50%. Cross sections were then obtained by 
applying bin-by-bin corrections to the inclusive jet distributions of the data sample in the 
variables nf et and rf^ x . The acceptance correction factors for the inclusive jet cross section 
da/dr] jet (T]^ x < 1.8) (a(r]^ x < r]^ ax )) were found to vary between 0.63 and 0.93 (0.59 and 
0.84). The dependence of these correction factors on the choice of parametrisations of the 
parton densities in the pomeron were found to be below ~ 20%, and are taken into account in 
the systematic uncertainty assigned to the measurements reported in the next section. 
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7 Results 



In this section, first the measured jet cross sections are presented and the uncertanties of the 
measurements discussed. These results are model- independent. Second, the expectations from 
non-diffractive processes are found not to account for the measurements. Third, the predictions 
from diffractive models are compared to the data and estimates of the momentum sum of the 
pomeron and of the relative contribution of quarks and gluons in the pomeron are extracted 
using solely the diffractive jet measurements. Fourth, the jet cross sections in photoproduction 
are combined with the measurements of the DIS diffractive structure function to constrain 
further the parton content of the pomeron. 



7.1 Jet Cross Sections 

The results for da j drf ei {j]^ x < 1.8) and a(r]^ x < ^ aa .) are presented in Figs. || and |3| and in 
Tables |I] and |2|. The differential cross section is flat as a function of if et . Since the measured 
jet cross sections refer to events with a large rapidity gap they include a contribution from 
double dissociation. The statistical errors of the measurements are indicated as the inner error 
bars in Figs. |2] and |3|. They are ~ 30% for da j drf et {j]^ x < 1.8) and constitute the dominant 
source of uncertainty. For a(r/^ x < ^ aa .) the statistical error increases from 8% to 20% as 
Vmax decreases. A detailed study of the systematic uncertainties of the measurements has been 
carried out |25|, [37|. The sources of uncertainty include the dependence on the choice of the 



parton densities in the pomeron, the simulation of the trigger, the cuts used to select the data, 
and the absolute energy scale of the cal jets . 



The following systematic uncertainties related to the rj max -c\xt were studied: the T]^ x variable 
in the data and the simulated events was recomputed after removing the CAL cells with rj > 3.25 
in order to check the dependence on the detailed simulation of the forward region of the detector, 
resulting in changes up to 13% for a(r]^ x < rf^^) and up to 18% for da j drf et {r\^ x < 1.8) 
(except at the most forward data point, where the statistics are small and the change amounts 
to 37%); the energy threshold in the computation of 77™^ for data and simulated events was 
decreased to 300 MeV, yielding changes up to 11% for da j drf et {r]^ x < 1.8) and up to 14% for 
a i r lmax < Vmax) (except at the most backward data point, where the statistics are small and 
the change amounts to 27%). 

The dominant source of systematic error is the absolute energy scale of the cal jets, known 
to within 5%, which results in a 20% error. The systematic uncertainties not associated with 
the energy scale of the jets were added in quadrature to the statistical errors and are shown 
as the total error bars. The additional uncertainty due to the energy scale of the jets is shown 
as a shaded band. The systematic uncertainties have large bin to bin correlations. They are 
to be understood as a conservative estimate of the error associated with each data point. An 
additional overall normalisation uncertainty of 3.3% from the luminosity determination is not 
included. 



7.2 Comparison to non-diffractive model predictions 

The contribution to the measured cross sections from non-diffractive processes was estimated 
using PYTHIA including resolved and direct processes. Had jets were selected in the generated^ 
events using the same jet algorithm as for the data and calculating 77^. as explained in section 

6 These generated events were analysed at the hadron level. 
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6. In PYTHIA the occurrence of a rapidity gap is exponentially suppressed and arises from a 
fluctuation in the pseudorapidity distribution of the final state hadrons. The calculations using 



MRSD_ |36] for the proton and GRV-HO |35|] for the photon parton densities are compared 
to the measurements in Figs. |] and 0. The non-diffractive contribution does not reproduce 
the measurements. For the measured dcr/dri^ et (r]^ x < 1.8) the non-diffractive contribution is 
close to the data only at the most forward measured point. For the remaining 7]i et range, the 
data are a factor between 2 and 7 above the expectations from non-diffractive processes. In the 
measured cr(i]^ x < ^ aa; ), the non-diffractive contribution as predicted by PYTHIA is smaller 
than the data by factors between 3 and 9. These comparisons, together with the features of 
the data shown in the previous section, demonstrate that the measured jet cross sections with 
a large rapidity gap cannot be accounted for by non-diffractive processes. However, in the 
discussion below, this non-diffractive contribution will be subtracted from the data. 



7.3 Comparison to diffractive model predictions 

The measured cross sections are compared to the predictions of the models for diffractive hard 
scattering mediated by pomeron exchange, as implemented in the POMPYT generator. The 
predictions have been obtained by selecting had jets in the generated events using the same jet 
algorithm as for the data and calculating r]^ x as explained in section 6. 

In a first step, the predictions of POMPYT, using the DL flux factor and the parametrisa- 
tions of the pomeron parton densities suggested theoretically (see section 3.1), and assuming 
Ejp = 1, are compared to the measured cross sections in Figs. § and |3|. For this initial com- 
parison, the contributions from non-diffractive and double dissociation processes have not been 
taken into account. As mentioned earlier (see section 3.1), the /i 2 dependence of the parton 
densities has been neglected and, hence, the argument of Ejp(/x 2 ) is omitted. The scale rel- 
evant for the measured jet cross sections is jj 2 ~ (E^) 2 . We start by discussing the results 
for da/dr] jet (r]^ x < 1.8). The shape of the predictions of POMPYT using a hard parton 
density compares well with the measured shape of the cross section. The shape predicted by 
POMPYT using a soft gluon density does not describe the data. The calculations based on 
a soft gluon density are smaller than the measurements by factors between 20 and 50. This 
type of parametrisation was already disfavoured by previous studies ||. The predictions using 
a hard quark density are too small by factors between 3 and 10, but those using a hard gluon 
density reproduce the measurements well. The predictions based on the IS flux factor lead to 
similar conclusions. For the integrated cross section o{j]^ x < rf^^), the measured shape is in 
each case described by the expectations from POMPYT, although the normalisation is incor- 
rect by a factor depending on the model. A soft gluon (hard quark) density yields a prediction 
which is too small by factors between 30 and 60 (5 and 10). A hard gluon density for the 
pomeron gives a good description of the data. Based on the samples of events of POMPYT, 
the measurements are sensitive to (3 values above approximately 0.3. Therefore, the data is not 
sensitive to a possible additional contribution due to a soft parton component in the pomeron. 
The data do not rule out a possible contribution from a super-hard parton component in the 
pomeron [15 — 17]. 

In principle, other processes could contribute to jet production with a large rapidity gap. 
For example, the proton may emit a tt + (instead of a pomeron), pi rif7c + , and a partonic 
constituent of the 7r + undergoes a hard interaction with the photon or its constituents. The 
contribution from this reaction to the data is expected to be small due to the power law 
decrease, ~ W~ A , for pion exchange. Monte Carlo calculations using POMPYT confirm these 
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expectations. 

In a second step, the data were compared to the predictions of POMPYT based on a pomeron 
consisting of both quarks and gluons but without assuming Sjp = 1. In addition, the con- 
tribution from non-diffractive processes and from double dissociation to the measured cross 
sections were taken into account. The non-diffractive contribution as predicted by PYTHIAQ 
was subtracted bin by bin from the data. The contribution from double dissociation for large- 
rapidity-gap events was estimated to be (15 ± 10)% []TU[. This contribution was assumed to be 
independent of 7]^ et and was also subtracted from the data. After the above subtractions, the 
data were compared with the predictions of POMPYT using the DL flux factor and allowing for 
a mixture of the hard gluon (6/?(l — /?)) and the hard quark (|/3(1— (3)) densities in the pomeron: 
a fraction c g for hard gluons and c q = 1 — c g for hard quarks. The overall normalisation of 
the POMPYT prediction was left as a free parameter: Sjp. For this study, the contribution 
to the cross sections and to Ejp from possible soft gluon and soft quark components has been 
neglected. 

For each value of c g , a one-parameter (£p) x 2_ fit to the measured da ' /drf et {rff^ < 1.8) was 
performed. The results are presented in Fig. [|. The thick solid line represents the value of Sjp 
for the minimum of the x 2_ fit for each value of c g and the shaded band represents the 1 a range 
around those minima. For c g = 1 (gluons only) the fit yields £jp = 0.5 ± 0.2 with xLin = 2-3 
for three degrees of freedom, while for c g = (quarks only) the fit yields Sjp = 2.5 ± 0.9 with 
Xmin = 2-8. The momentum sum rule (Sjp = 1) is approximately satisfied for 0.2 < c g < 0.6 
(statistical errors only). Note that for this estimate the DL form for the pomeron flux factor 
was assumed. 

This comparison of cross sections for jet production with a large rapidity gap between data 
and model predictions is subject to the following uncertainties: 

• The jet cross sections obtained from the Monte Carlo calculations presented here are 
leading order calculations. In these calculations, a s {^ 2 ) and the parton densities in the 
proton and the photon are evaluated at /i 2 = f>%. These computations may be affected by 
higher order QCD corrections, which are expected to change mainly the normalisation (K- 
factor). The agreement between the PYTHIA calculations of the inclusive jet differential 



cross sections and the measurements p5| indicate that in the case of the non-diffractive 
contribution the i-T-factor is close to 1, within an uncertainty of ~ 20%. The fT-factor in 
the case of POMPYT is expected to be similar (with a similar uncertainty), as the same 
hard subprocesses are involved in the calculation of jet cross sections. 

The amount of the non-diffractive contribution to the measured cross section was modelled 
using PYTHIA with some choices for the proton and photon parton densities. This 
contribution is more sensitive to the choice of photon parton densities. 

The uncertainty in the estimation of the contribution from double dissociation. 

The POMPYT model for diffractive hard scattering assumes factorisation of the hard 
process with respect to the soft diffractive reaction. The extent to which this assumption is 
valid has to be determined experimentally through a detailed comparison of measurements 
for different reactions (see next section). 



7 These calculations give a good description of the inclusive jet differential cross sections (without the large- 
rapidity-gap requirement) in the range — 1 < rf et < 1 ]25j. 
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• The pomeron flux factors adopted in the various models are based on different assumptions 
for the t and xp dependences which are obtained from data on soft diffractive hadronic 
processes. The uncertainty in the procedure used to extract the flux is about 30%. 

The differences between the results obtained in each of the studies listed above and the 
central values were combined in quadrature to yield the theoretical systematic uncertainties 
(not shown in Fig. f|) of the fitted values of Sjp. These uncertainties were then added in 
quadrature with the statistical and systematic uncertainties of the measurements resulting in 
the following ranges at the 1 a level: 1.4 < Tip < 3.8 for c g = and 0.3 < Tp < 0.9 for c g = 1. 
The range in Tp assumes the DL convention for the pomeron flux factor. This normalisation 



has recently been discussed by Landshoff ]38] who concludes that the normalisation is arbitrary 
up to a multiplicative factor A. If the normalisation is changed by a factor A, the range of the 
momentum sum is given by 1A/A < Tp < 3.8/ A for c g = and by 0.3/ A < Tp < 0.9/ A for 
c g = 1. 

In summary, the comparison of model predictions with the jet cross section measurements 
favours those models where the partonic content of the pomeron has a hard contribution. Given 
the uncertainties mentioned above and the DL convention for the normalisation of the pomeron 
flux factor, the data can be reproduced by a pomeron whose partonic content varies between a 
pure hard quark density with momentum sum given by 1.4 < Tp < 3.8 and a pure hard gluon 
density with 0.3 <Tp < 0.9. 

7.4 The gluon content of the pomeron 

The HERA experiments have recently presented the first measurements of the diffractive struc- 
ture function in DIS || [T(J. The results show that the quark densities in the pomeron have a 



hard and a soft contribution. Assuming the DL form for the pomeron flux factor, the DIS data 
do not favour a pomeron structure function which simultaneously fulfils the momentum sum 
rule and consists exclusively of quarks. 

If the pomeron parton densities are universal and describe both DIS and photoproduction 
processes, the DIS results together with the photoproduction data further constrain the partonic 
content of the pomeron. The measured diffractive structure function in DIS (F^ 3 \f3,Q 2 ,xp)) 



TOfl can be used to extract the contribution of the quarks to the momentum sum (Tp q (Q )). 



The integral of over xp and j3 is proportional to Tp q (Q 2 ) 

IP max /*1 r~>fQ\ n r% 

/ dxjp / d(3 F 2 Di3 \(3,Q 2 , XlP ) = k f ■ Tp q (Q 2 ) ■ I flux (9) 

where Ifi ux is the integral of the pomeron flux factor over t and over the same region in 
xjp, and kf is a number which depends on the number of flavours assumed (5/18 for two 
flavours and 2/9 for three flavours). For the left-hand side of Eq. (9), the parametrisation 
of F 2 D{3) ((3,Q 2 , Xjp) obtained in |I| was used. The integral was performed over the range 
6.3 • 10 -4 < xp < 10~ 2 of the ZEUS DIS measurements. The DL form for the pomeron flux 
factor was used to compute Ifi ux for the right hand side of Eq. (9). This procedure yields an 
estimate of T,p q (Q 2 ): 0.32 ± 0.05 (0.40 ± 0.07) for two (three) flavours. These estimates are 
based on a parametrisation of F^ 3 \/3,Q 2 ,xp) which was determined in the large (3 region 
{0.1 < (3 < 0.8) and is assumed to be valid for the entire region < (3 < 1. In the range of 
Q 2 where the DIS measurements were done, 8 GeV 2 < Q 2 < 100 GeV 2 , the pomeron structure 
function is approximately independent of Q 2 and, thus, the estimated Sjp g (Q 2 ) does not depend 
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upon Q 2 . It should be noted that the scales at which the parton densities in the pomeron are 
probed in DIS, Q 2 , and in photoproduction, /i 2 , are comparable. The estimate from DIS 
imposes a constraint on the Sjp — c g plane which, combined with the estimates obtained in the 
preceding section, restricts the allowed ranges for E#> and the relative contributions of quarks 
and gluons (c g ). The DIS constraint, which can be written as Sjp ■ (1 — c g ) = 0.32 (0.40) for 
the two choices of the number of flavours, is included in Fig. ^ (the dark shaded area represents 
the uncertainty in this constraint). Combining the estimates from photoproduction (thick solid 
line) and DIS yields 0.5 < Sjp < 1.1 and 0.35 < c g < 0.7 (statistical errors only). 

These results are subject to the uncertainties listed at the end of section 7.3. The allowed 
range for Sjp which results from the combination of the DIS and photoproduction measure- 
ments was evaluated for each source of systematic uncertainty. Taking into account all the 
uncertainties mentioned, the comparison between the DIS and photoproduction measurements 
gives 0.4 < Sjp < 1.6 for the momentum sum of the pomeron assuming the DL convention for 
the flux. If the normalisation of the pomeron flux factor is changed by a multiplicative factor 
A, the allowed range of the momentum sum is given by OA/A < Sjp < 1.6/ A. 

It should be noted that the evaluation of the c g range allowed by the DIS and photoproduction 
measurements is not affected by the normalisation of the pomeron flux factor or the uncertainty 
on the double dissociation contribution since they cancel out in the comparison^ Taking into 
account the remaining uncertainties the combination of the DIS and photoproduction data 
gives 0.3 < c g < 0.8. This result does not depend on the validity of the momentum sum rule 
for the pomeron. 

8 Summary and conclusions 

Measurements of ep cross sections for inclusive jet photoproduction with a large rapidity gap in 
ep collisions at yfs = 296 GeV using data collected by the ZEUS experiment in 1993 have been 
presented. The measured jet cross sections are compared to perturbative QCD calculations of 
diffractive hard processes and allow a model dependent determination of the parton content of 
the pomeron. The measurements require a contribution from a hard momentum density of the 
partons in the pomeron. This result is consistent with the observations of the UA8 Collaboration 
made in pp collisions. When the measured jet cross sections are combined with the results on the 
diffractive structure function in deep inelastic scattering at HERA, first experimental evidence 
for the gluon content of the pomeron is found. This evidence is independent of the normalisation 
of the flux of pomerons from the proton and does not rely on assumptions on the momentum 
sum of the pomeron. The data indicate that between 30% and 80% of the momentum of the 
pomeron carried by partons is due to hard gluons. 
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da/dr]i et (ri^ x < 1.8)±stat.±syst.±syst. £f et -scale 


range 


fin nbl 

V 111 F 1 -"/ 


(-1,-0.5) 


44 ± 12 ± 8 ±1° 


(-0.5,0) 


39 ±11 ±5 If 


(0,0.5) 


37 ± 10 ± 10 tjj 


(0.5,1) 


16 d= 6 ± 7 ±| 



Table 1: Measured differential ep cross section da /drf et {r^ x < 1.8) for inclusive jet production 
for E J T et > 8 GeV and in the kinematic region Q 2 < 4 GeV 2 and 0.2 < y < 0.85. The r] jet 
ranges used for the measurements are shown. The cross sections are given at the centre of each 
rf et bin. The statistical and systematic errors are also indicated. The systematic uncertainties 
associated with the energy scale of the jets are quoted separately. The overall normalisation 
uncertainty of 3.3% is not included. 



„0 
Imax 


cr(r)max < r^ a J±stat.±syst.±syst. E J T et -scale 
(in pb) 


1.0 


27.3 ±5.5 ±9.9 t|| 


1.5 


51.7 ±8.5 ±9.2 ±^i 7 


1.8 


67 ± 10 ± 12 t\l 


2.0 


98 ± 12 ± 14 t™ 


2.2 


145 ± 15 ± 29 til 


2.4 


208 ± 18 ± 43 tfl 



Table 2: Measured integrated ep cross section <r{r]^ x < ^ a;c ) for inclusive jet production for 
E J T et > 8 GeV and -1 < r] jet < 1 in the kinematic region Q 2 < 4 GeV 2 and 0.2 < y < 0.85. The 
statistical and systematic errors are also indicated. The systematic uncertainties associated 
with the energy scale of the jets are quoted separately. The overall normalisation uncertainty 
of 3.3% is not included. 
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Figure 1: (a) The scatter plot of Mjf z versus rj^ x for the sample of events with at least 
one cal jet fulfilling the conditions E^l al > 6 GeV and — 1 < rfjjj < 1; (b) the distribution 
of rj™ l ax for the events with Mjf ' < 30 GeV along with the predictions of PYTHIA (shaded 
area) and POMPYT with a hard gluon density in the pomeron (solid line). The predictions 
are normalised to the number of data events above and below rf^ ax = 2.5, respectively; (c) the 
distribution in M™ 1 for the events with rf^ ax < 1.8 together with the prediction of POMPYT 
with a hard gluon density in the pomeron (solid line) normalised to the number of data events; 
(d) the distribution in W cal for the events with rf^ ax < 1.8 and the prediction of POMPYT as 
in (c). 
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Figure 2: Measured differential ep cross section da/drf et (r]!^x < 1-8) for inclusive jet produc- 
tion for E?p et > 8 GeV in the kinematic region Q 2 < 4 GeV 2 and 0.2 < y < 0.85 (dots). The 
measurements are not corrected for the contributions from non-diffractive processes and double 
dissociation. The inner error bars represent the statistical errors of the data, and the total 
error bars show the statistical and systematic errors —not associated with the energy scale of 
the jets— added in quadrature. The shaded band displays the uncertainty due to the energy 
scale of the jets. For comparison, POMPYT predictions for single diffractive jet production 
(e + p — > e + p + jet + X r ) using the DL flux factor for direct plus resolved processes for various 
parametrisations of the pomeron parton densities (hard gluon, upper solid line; hard quark, 
middle solid line; soft gluon, lower solid line) are also shown. The GS-HO photon parton densi- 
ties have been used in POMPYT. The contribution from non-diffractive processes is exemplified 
by the PYTHIA predictions using MRSD_ (GRV-HO) for the proton (photon) parton densities 
(dashed line). 22 
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Figure 3: Measured integrated ep cross section o{j]^ x < fj^ax) f° r inclusive jet production for 
-1 < 7] jet < 1 and E 3 T et > 8 GeV in the kinematic region Q 2 < 4 GeV 2 and 0.2 < y < 0.85 
(dots). The measurements are not corrected for the contributions from non-diffractive processes 
and double dissociation. The inner error bars represent the statistical errors of the data, and 
the total error bars show the statistical and systematic errors —not associated with the energy 
scale of the jets— added in quadrature. The shaded band displays the uncertainty due to the 
energy scale of the jets. For comparison, PYTHIA and POMPYT calculations (for the same 
conditions as in Fig. 2) are included. 
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Figure 4: The plane of the variables £jp (momentum sum) and c g (relative contribution of 
hard gluons in the pomeron). The thick solid line displays the minimum for each value of c g 
obtained from the \ 2 fit (the shaded area represents the 1 a band around these minima) to 
the measured da/drf et (r]^ x < 1.8) using the predictions of POMPYT. The constraint imposed 
in the Sjp — c g plane by the measurement of the diffractive structure function in DIS (F®^) 
[10] for two choices of the number of flavours (upper dot-dashed line for Y>jp q = 0.40 and lower 
dot-dashed line for Ejp 9 = 0.32) is also shown. The horizontal dashed line displays the relation 
SiP = 1. 
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